Universality in solar flare and earthquake occurrence 
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Earthquakes and solar flares are phenomena involving huge and rapid releases of energy char- 
acterized by complex temporal occurrence. By analysing available experimental catalogs, we show 
that the stochastic processes underlying these apparently different phenomena have universal prop- 
erties. Namely both problems exhibit the same distributions of sizes, inter-occurrence times and 
the same temporal clustering: we find afterfiare sequences with power law temporal correlations as 
the Omori law for seismic sequences. The observed universality suggests a common approach to the 
interpretation of both phenomena in terms of the same driving physical mechanism. 
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Solar flares are highly energetic explosions |lj from ac- 
tive regions of the Sun in the form of electromagnetic ra- 
diation, particle and plasma flows powered by strong and 
twisted magnetic fields. Since they cause disturbances on 
radio-signals, satellites and electric-power on the Earth, 
much interest has been devoted in the last years to space 
weather forecasts [2J, y] • Recent studies have shown that 
solar flares also affects the Sun's interior, generating seis- 
mic waves similar to earthquakes 4j. Actually, despite 
having different causes, solar flares are similar to earth- 
quakes in many respects, for example in the impulsive 
localised release of energy and momentum and their huge 
fluctuations The analogy with earthquake occur- 

rence is also supported by the observation of power laws 
0, 0, in the distribution of flare sizes, P(s), related 
to the Gutenberg-Richter law for the earthquake mag- 
nitude distribution. Various interpretations have been 
proposed for these power law distributions ranging from 
Magneto- Hydro-Dynamics mpM to turbulence |lfj up to 
Self Organized Criticality (nUllll^. A better under- 
standing of solar flares and coronal mass ejections from 
the Sun requires knowledge of the structural details of 
these events and their occurrence in time. This could 
greatly improve the prediction of violent space weather 
and the understanding of the physical processes behind 
solar events. 

Here we present evidence that the same empirical laws 
widely accepted in seismology characterize, surprisingly, 
also the size and time occurrence of solar flares. In par- 
ticular the same temporal clustering holds both for earth- 
quake, where it is known as the Omori law, and solar flare 
catalogues: a mainflare triggers a sequence of afterflares. 
The evidence of a universal statistical behaviour suggests 
the possibility of a common approach to long term fore- 
casting and rises as well deep questions concerning the 
nature of the common basic mechanism. 



A statistical approach to earthquake catalogues has 
revealed a scale invariant feature of the phenomenon, as 
indicated by power law distributions for relevant phys- 
ical observables [Til llij . such as the seismic moment 
distribution of earthquakes, P{s) ~ s~", where the ex- 
ponent a G [1.6, 1.71 is essentially the same in different 
areas of the world [l6|. This relation corresponds to the 
Gutenberg-Richter law for the distribution of the earth- 
quake magnitude M via the relation M = 2/3 log(s) — K, 
where if is a constant ^3 ■ R is a l so observed that earth- 
quakes tend to occur in clusters temporally located after 
large events: the Omori law states that the number of 
aftershocks at time t after a main event, NaU ), decays 
as a power law Na(£) ~ t~ v with p ~ 1 [18J. Finally, 
the distribution of intertimes between consecutive earth- 
quakes, P(At), is not a simple power law, but has a non 
trivial functional form which, like the other quantities 
mentioned before, is essentially independent of the geo- 
graphical region or the magnitude range considered |l9| . 
These observations suggest that P(At), NA(t) and P(s) 
are distinctive features of earthquakes and, thus, funda- 
mental quantities for a probabilistic analysis of the phe- 
nomenon characterizing its amplitude and time scales. 

In this letter we analyse several catalogues of solar 
flares and compare them with the Southern California 
catalogue for earthquakes. Since emissions at differ- 
ent wavelengths are related to different radiation mech- 
anisms, we present a comparison among solar data from 
X-ray observations in three different energy ranges and 
different periods of solar cycle, by using on-line available 
flare catalogues. More specifically, Soft X-ray data in the 
(1.5-2.4) keV and (3.1-24.8) keV ranges are obtained from 
the Geostationary Operational Environmental Satellite 
(GOES) systems HJ. We consider N e = 21567 events 
from January 1992 to December 2002 covering the entire 
11-years solar cycle. Solar flares in the Hard X-ray range 
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FIG. 1: (Color online) The number distribution, n(At), of 
intertimes, At, between consecutive events in solar flare and 
earthquake catalogues. Solar data refer to X-ray observations 
in three different energy ranges covering different periods of 
solar cycle: soft X-ray data in the (1.5-2.4) keV and (3.1-24.8) 
keV ranges from the GOES catalogue (■ symbols); hard X- 
ray (>25 keV) from the BATSE catalogue (O symbols); in- 
termediate X-ray (10-30 keV) from the WATCH catalogue (□ 
symbols) . Earthquake intertimes data are from the California 
catalogue for events with magnitude M > 2 (• symbols). 



(>25 keV) are obtained from the Burst and Transient 
Source Experiment (BATSE) that gives N e — 6658 flares 
from April 1991 to May 2000 Finally N e = 1551 

events from January 1990 to July 1992 in the Interme- 
diate X-ray (10-30 keV) range are analysed from the 
WATCH experiment [22(. Many earthquake catalogues 
exist and the universality of their statistical features has 
already been established Thus, for clarity, we 

only consider here the Southern California earthquake 
catalogue [23| which has N e = 88470 events with magni- 
tude M > 2 in the years from 1967 to 2002. 

The intertime distribution has been already investi- 
gated both for earthquakes 0] and solar flares [24] . The 
intertime, At, is the time between the start of a flare 
(or an earthquake) and the start of the next one as re- 
ported in the above catalogues. Here, for a catalogue 
with N e data, we count the number of events n(At) 
having intertime between At and At + X/N e , where A 
is a constant setting the binning of raw data. This is 
the statistically relevant quantity to consider [25}, since 
n(At)/\ — > P(At) in the limit N e 00, and, thus, in 
the following we refer to n(At). Here we choose A such 
that \/N e = Ihsec for the California catalogue and use 
the same value for all catalogues. Fig.l shows the in- 
tertime distributions, n(At), for the three different solar 



FIG. 2: (Color online) The number distribution, n(At), 
of intertimes, At, between consecutive solar flares for the 
GOES catalog. Data from GOES8 and GOES 10 satel- 
lites correspond to the minimum of the solar cycle (from 
9/1/1995 to 12/31/1996) and the maximum (from 1/1/2000 
to 12/31/2003) respectively. 



flares data sets and for the California earthquake cata- 
logue. Solar flares data scale one on top of the other to a 
very good approximation and, interestingly, they all ap- 
pear to collapse, within statistical errors, on the same non 
trivial distribution function of earthquake intertimes. In 
particular, this data collapse is obtained without rescal- 
ing At by any suitable factor: the intertime duration, 
At, is expressed in the same units (seconds) for all data 
sets. Thus, Fig.l shows that the same intertime distribu- 
tion function and, surprisingly, even the same time range 
characterize these apparently different physical processes 
in the magnitude range of the above catalogues. 

The scaling behaviour of n (At) for different solar 
phases is a widely debated sub ject |2rj I27I Ejf and a 
dependence on the solar phase [26J has been observed 
also in the case of Coronal Mass Ejection [27J. The re- 
sult for flares has been obtained by taking into account 
only events with a peak flux larger than 1.4 x 10~ 6 Wm~ 2 
(class CI). We have then considered separately data from 
the GOES catalog corresponding to maximum and mini- 
mum solar activity and used the same binning procedure 
as Fig.l. In order to take into account the different level 
of background X-ray flux, we have set different thresh- 
olds for different phases: events greater of class CI in the 
maximum and class Bl (10 -7 x Wm~ 2 ) in the minimum 
phase 29]. Fig. 2 shows that data from different phases 
fall on the same universal curve. 

The other crucial quantity to be investigated is the dis- 
tribution of flares sizes, P{s), i.e., the distribution of the 
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FIG. 3: (Color online) The distribution, n(s/so), of flare 
peak intensity, from the same catalogues of Fig.l, and of seis- 
mic moments from the California earthquake catalogue. We 
set A/JV e = 1 for the California catalogue. As for the inter- 
times of Fig.l, comparison between the size distributions of 
earthquakes and solar flares shows very good agreement. The 
universal distribution is well fitted by a power law with ex- 
ponent a — 1.65 ± 0.1 (shown as a solid line in the picture). 
s = 10" 7 W/m 2 for GOES, s = 600 cmnts/(sec 2000cm 2 ) 



for BATSE, s = 3000 cnt/bin for WATCH, s = 30 
Nm for earthquakes. 
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flare peak intensity, s, from the above catalogues. This 
is compared with the earthquake seismic moments dis- 
tribution, P(s), from the Southern California catalogue. 
In order to normalize the different units and experimen- 
tal ranges used in each catalogue, here we scale the val- 
ues, s, of each data set by a given constant amount So- 
Then, we calculate the number of events n(s/so) with 
sizes between s/sq and s/sq + X/N e . Excellent data col- 
lapse is observed in Fig. 3 with all data fitted over almost 
three decades by a power law n(s/so) ~ (s/so) _a with 
an exponent a = 1.65 ± 0.1, in agreement with previ- 
ous results on solar flares 0, 0, and earthquakes 0] . 
Therefore, in analogy to earthquakes, from the above ob- 
servations it is possible to introduce a Richter scale for 
flares where their "magnitude" , M, is defined via the re- 
lation: M(s) = 2/31og(s) ~ Kp, where Kp is a constant. 
In terms of the magnitude the data from flares catalogs 
are therefore found to follow the Gutenberg Richter law 
introduced for earthquakes. 

Further evidence of structural similarities in the statis- 
tics of the two phenomena is given by the analysis of cor- 
relations between events within each of these catalogues. 
It would be interesting to compare the time correlation 
between main-events and the sequence of their after- 



events, as in the Omori law. We define a "main-event" as 
an event with magnitude M > M ma i n ; its "after-events" 
are the following events with M cut < M < M ma i n , where 
M cut is a cutoff for small background events. The basic 
difference with the usual definition used in seismology 
is that an event with M < M ma i n considered as "after- 
shock" may instead be an independent event totally un- 
related to the preceding "mainshock" . Furthermore, an 
event with M > M ma i n considered as "mainshock" may 
have been triggered by a previous larger event. Despite 
these differences, our definition can be straightforwardly 
applied to flare catalogues too and tends to the stan- 
dard one for large enough M ma i n and M cut : here we fix 
M cu t = M ma i n — 2.5. In Fig. 4 we show the number of 
"after-event", n^i(t), at time t after a "main-event", for 
all the mentioned data sources. Interestingly, the time 
correlation function, riyi(t), has the same functional form 
in all catalogues. A power law nA(t) ~ 1/t (straight line 
in the picture), as the Omori law, fits the data. The 
results are quite robust with respect to changes of the 
parameter M ma i ni provided that M ma i n is large enough 
as previously explained. We apply the same procedure 
to analyse the rate of occurrence of events leading up to 
a main event and observe that also " fore-flares" follow 
the same power law behaviour (Omori law) as foreshocks 
|30j . even if more symmetrical behaviour is observed in 
the flare case. 

Fig.s 1, 3 and 4 indicate that the statistical properties 
of size and time scales of solar flares (independently of 
the energy range and the temporal location in the so- 
lar cycle of the X-ray radiation) and earthquakes are 
essentially the same within current statistical accuracy. 
It is tempting to look at the observed universality in 
the perspective of the theory of critical phenomena. In 
the past analogy between the two phenomena was pro- 
posed on the basis of the same theoretical model \?>\\ . 
Here we follow a completely different approach: we di- 
rectly compare experimental catalogs, we observe uni- 
versal behaviour and therefore we propose the presence 
of a common driving physical mechanism. Most earth- 
quakes occur where the elastic energy builds up owing 
to relative motion of tectonic plates. Schematically, as 
the friction locks the sliding margins of the plates, en- 
ergy load increases. When elastic stress overcomes the 
threshold of frictional resistance, it is relaxed causing 
the occurrence of an earthquake. This "stick-slip" be- 
haviour redistributes the stress-energy field in the crust 
generating new earthquakes where and when the local 
slipping threshold is exceeded. A quantitative predic- 
tion on the aftershocks decay cannot be derived by sim- 
ple stress transfer but can be obtained in terms of a 
state variable constitutive formulation, where the rate 
of earthquakes results from the applied stressing history 
|32j |. This formulation gives account for long-range cor- 
relations between earthquakes affecting the shape of the 
whole intertime distribution. 
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tablished methods of seismic forecasting j.3 l| . 
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FIG. 4: (Color online) The correlation function, nA(i), i.e., 
the number of "after-event" at time t after a "main-event" 
for the same catalogues of Fig.l. To find the best collapse, 
data from different catalogues are rescaled by a given amount 
t (t = 700 sec for GOES, t = 60 sec for BATSE, to = 1 
sec WATCH and to = 43 sec for California earthquakes). As 
for Fig.l and Fig.3, the "aftershock" rate of occurrence for 
earthquakes and solar flares scale very well. For comparison, 
we also plot an Omori power law n,A(t/to) ~ to/t (solid line). 



The universal scaling of Figures 1-4 suggests a sim- 
ilar physical mechanism at the basis of solar flares oc- 
currence. Flares and X-ray jets arise in active solar re- 
gions where magnetic flux emerges from the solar interior 
and interacts with ambient magnetic field. These inter- 
actions are thought to occur in electric current sheets 
separating regions of opposite magnetic polarity. The 
dynamics and energetics of these sheets are governed by 
a complex magnetic field structure . Opposite fluxes 
lead to rearrangement of field lines building up magnetic 
stress up to a breaking point where magnetic energy is 
released in a flare via magnetic reconnections. The ob- 
served temporal clustering of Fig. 4 shows that the rate 
of flare occurrence decreases in time as a power law af- 
ter a " mainflare" . Since the same behaviour is found 
for seismic sequences, here we propose that the mech- 
anism at the basis of seismic energy redistribution can 
be responsible for "afterflare" occurrence. In particular 
magnetic stress transfer in Solar Corona plays the role 
of elastic stress redistribution on the Earth crust. As a 
consequence the state-rate formulation can be gener- 
alized to solar flares, namely the flare triggering depends 
on the entire history of magnetic stresses. Beyond issues 
of fundamental science, the present results can also have 
very practical consequences as, for instance, to improve 
the prediction of violent space weather by applying es- 
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